Nuclear Spin-Lattice Relaxation in One-Dimensional Heisenberg Ferrimagnets: 

Three-Magnon versus Raman Processes 
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Nuclear spin-lattice relaxation in one-dimensional Heisenberg ferrimagnets is studied by means of 
a modified spin-wave theory. We consider the second-order process, where a nuclear spin flip induces 
virtual spin waves which are then scattered thermally via the four-magnon exchange interaction, as 
well as the first-order process, where a nuclear spin directly interacts with spin waves via the hy- 
perfine interaction. We point out a possibility of the three-magnon relaxation process predominating 
over the Raman one and suggest model experiments. 



PACS numbers: 75.10.Jm, 75.50.Gg, 76.50.+g 

Design of molecule-based ferromagnets is a challenging 
topic in materials science [1] and numerous quasi-one- 
dimensional ferrimagnets have been synthesized in this 
context. Bimetallic chain compounds [2,3] are early ex- 
amples, one of which has indeed accomplished the three- 
dimensional ferromagnetic order. Another approach [4] 
consists of bringing into interaction metal ions and stable 
organic radicals. Genuine organic ferrimagnets [5,6] were 
also synthesized. The tetrameric bond-alternating chain 
compound Cu(C5H4NCl)2(N3)2 [7] and the trimeric in- 
tertwining double-chain material Ca3Cu3(PC>4)4 [8] are 
distinct ferrimagnets of topological origin. 

Theoretical understanding of one-dimensional ferri- 
magnets has also significantly been developed [9-14], but 
very little is known about their dynamic properties yet 
[15,16]. Although the one-dimensional spin dynamics is 
a modern topic of great interest [17], systematic calcula- 
tion of the nuclear magnetic relaxation is still absent. In 
the late sixties, Pincus and Beeman [18] formulated the 
nuclear spin-lattice relaxation rate for Heisengerg ferro- 
magnets and antiferromagnets by means of the spin- wave 
theory and pointed out a significant contribution of the 
three-magnon scattering to the relaxation rate. However, 
their theory is not effective in one dimension but is valid 
at the onset of the three-dimensional long-range order. 
In such circumstances, modifying the conventional spin- 
wave theory, we make our first attempt at systematically 
describing the one-dimensional nuclear spin dynamics on 
the basis of the spin-wave picture. We predict that the 
three-magnon relaxation process should predominate over 
the Raman one at high temperatures and weak fields in 
one- dimensional Heisenberg ferrimagnets. 

First of all our scheme [12,19] of modifying the spin- 
wave theory is distinct from the original idea proposed 
by Takahashi [20] and Hirsch et al. [21]. In the original 
way of suppressing the divergence of the sublattice mag- 
netizations, a Lagrange multiplier is first introduced and 
then the effective Hamiltonian is diagonalized subject to 
zero staggered magnetization. The thus-obtained energy 
spectrum depends on temperature and fails to describe 
the Schottky peak of the specific heat [22]. In order to 
obtain better thermodynamics, we first diagonalize the 



Hamiltonian, keeping the dispersion relations free from 
temperature, and then introduce a Lagrange multiplier 
so as to minimize the free energy. This scheme works well 
for ferrimagnets, as demonstrated in Fig. 1. Both the ap- 
proaches well reproduce the magnetic susceptibilty, but 
our new shceme is much better than the original one at 
describing the specific heat. The nonvanishing specific 
heat at high temperatures, that is, the endlessly increas- 
ing energy with temperature, is a fatal weak point of the 
original scheme. 

We consider Heisenberg ferrimagnetic chains of alter- 
nating spins S and s, as described by the Hamiltonian 
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FIG. 1. Modified spin-wave calculations of the specific heat 
and the magnetic susceptibility as functions of temperature 
for the spin-(|, ~) ferrimagnetic Heisenberg chain. The orig- 
inal (Takahashi) and our new (Yamamoto) schemes are com- 
pared with numerical findings (Exact) [10]. 
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FIG. 2. Illustration of the elementary nuclear spin-lattice 
relaxation processes. Solid arrows designate ferromagnetic or 
antiferromagnetic spin waves inducing a nuclear spin flip (x). 
Broken arrows denote the four-magnon exchange interaction, 
(a) The direct process, (b) The first-order Raman process, (c) 
The first-order three-magnon process, (d) The second-order 
three-magnon process, where q = — k,A = ks — — fci. 
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Introducing the Fourier-transformed Holstein-Primakoff 
magnon operators a k and b k for each sublattice, we 
expand the Hamiltonian with respect to l/S as H = 
H 2 +n 1 +Ho+0(S l ), where we assume that O(S) = O(s) 
and Hi contains the 0(S l ) terms. The 0(S^ 1 ) contribu- 
tions are neglected in the following. The hyperfine inter- 
action is assumed to be of the dipolar type as 

W h f = 5MbS-7n^ + 



&n + B n s n ) + A z n S^ + B z 



(2) 



where A° and B° are the hyperfine coupling tensors be- 
tween the nuclear and electronic spins. 

Since Hq and Hhf are both much smaller than Hi, we 
may regard Ho + Hm = V as perturbative interactions to 
the linear spin-wave system 

H 2 +H x = -2SsJN -(S + s)JN - gn B H(S - s)N 

(3) 



where at 



a k cosh6 k + b k smh9 k and f3 k = afcsinh^fe + 

b\cosh8k, provided tanh20fc = 2y/Ss cos(fc/2)/ (S + s), 
create the ferromagnetic and antiferromagnetic spin 
waves of momentum fc, respectively, whose energies are 
given by u> k = u> k + <r(S — s) — og\i B H with u> k = 
[{S - s) 2 + 4Sssin 2 (fc/2)] 1 / 2 . We calculate the thermal 
distribution functions, (a k a k ) = n k and (/3t/3fc) = nt, 
by enforcing the zero staggered-magnetization constraint 
[12]. If we consider up to the second-order perturbation 
with respect to V, the probability of a nuclear spin being 
scattered from the state of I z = m to that of I z = m + 1 
is expressed as 
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m-E f ), 
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time is then given by T\ = (I— m)(I+m+l)/2W. Equa- 
tion (4) contains various elementary relaxation processes, 
which arc illustrated in Fig. 2. Due to the electronic- 
nuclear energy conservation, the direct process, involving 
a single spin wave, is rarely of significance. Within the 
first-order mechanism, where a nuclear spin directly in- 
teracts with spin waves via the hyperfine interaction, the 
three-magnon relaxation rate is much smaller than the 
Raman one [18]. However, the first-order relaxation rates 
are generally enhanced through the second-order mecha- 
nism, where a nuclear spin flip induces virtual spin waves 
which arc then scattered thermally via the four-magnon 
exchange interaction. We consider the leading second- 
order relaxation, that is, the exchange-scattering-induced 
three-magnon process, as well as the first-order relax- 
ation. The second-order relaxation of the Raman type 
originates in two virtual spin waves within the present 
model and is therefore negligible in our argument. 

Because of the significant difference between the elec- 
tronic and nuclear energy scales, hui^ -C JojZ < Ju k , 
the intraband spin-wave scatterings determine the Ra- 

(2) 

man relaxation rate 1/T| , whereas both the intra- 
band and interband spin-wave scatterings are relevant 
to the three-magnon relaxation rate 1/Tf . Assuming 
the Fourier components of the coupling constants to have 



little momentum dependence as J2 n e A n n 
and£ r 



^ikn 



we obtain 
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where i and / designate the initial and final states of the 
electronic-nuclear spin system described by the unper- 
turbed Hamiltonian. The nuclear spin-lattice relaxation 



where k 2 and fcf are given by u> ki 
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sinh6*fe, and the four-magnon exchange interaction Vi is given as 
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We calculate the case of (S,s) = (|, |), which is rel- 
evant to several major materials [2-4], assuming that 
A~/2 ~ A z = A and T^/2 ~ B z = B. Figure 3 shows 
1/Ti as a function of temperature and an applied mag- 
netic field. All in all, the exchange-scattering-enhanced 
three-magnon relaxation rate l/T^ grows with increas- 
ing temperature and decreasing field and ends in the lead- 
ing contribution to 1/Ti. The temperature dependence 
of 1 /Ti is determined by unless temperature is suffi- 
ciently high. As temperature increases, is decreasing 
at k ~ but otherwise increasing [23]. \duj^ / dk\~_ k _ 

makes the fci ~ components predominant in the mo- 
mentum summation (5), while such predominance arising 
from {dui^ /dk\~^ k _ is suppressed in the double summa- 

(2) 3 (3) 

tion (6) . 1 /Ti and 1 /T{ are hence decreasing and in- 
creasing, respectively, with increasing temperature and 
therefore, the Raman-to-three-magnon crossover may be 
detected at low fields. 



The field dependence of 1 /Ti is also useful in detecting 
the crossover. At moderately low temperatures and weak 
fields, fruj-N <C fceT "C J, we may evaluate Eq. (5) as 

1 _ 2( WB h lN ) 2 (AS-Bs) 2 

-gn&H/k B T ry- ( \ 

jf)~ irhSs(S-s)J °V2fc B T/' 

(10) 

where Kq is the modified Bcssel function of the sec- 
ond kind and behaves as K (hLo N /2k B T) ~ 0.80908 - 
ln(fi,WN/fcBT) for K,wn <C k B T. Considering the relation 

(2) 

= TN-ffj w e find that the field dependence of 1/T{ 
is logarithmic at weak fields and turns exponential with 
increasing field. Since 1/Ti exhibits much stronger ini- 
tial field dependence, the three-magnon relaxation can be 
detected at weak fields. We plot the crossover points in 
Fig. 4. The three-magnon relaxation process generally 
predominates over the Raman one at high temperatures 
and weak fields. Since the ferrimagnctic nuclear spin- 
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FIG. 3. Modified spin-wave calculations of the tempera- 
ture (the left three) and field (the right three) dependences of 
the nuclear spin-lattice relaxation rate at A/B — 0.1, where 
\/T[ 2 ^ and 1/T^ are plotted by dotted and broken lines, re- 
spectively, while 1/T 1 (2) + 1/t[ 3) = 1/Tj, which is observable, 
by solid lines. 
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FIG. 4. The crossover point, where 1/t[ 2) = 1/T^ 3) , as a 
function of temperature and an applied field. 

lattice relaxation is very sensitive to another adjustable 
parameter A/B, that is, the location of the probe nuclei, 
any T\ measurement should be accompanied by detailed 
structural analyses. At the special location of A/B ~ 
(dg/dsY — s/S, where ds (d s ) is the distance between 
the nuclear and larger (smaller) electronic spins, the lead- 
ing ferromagnetic spin waves are almost invisible to the 
nuclear spins [16,23] and therefore 1/Ti is too small to be 
measured. The three-magnon-dominant relaxation rate 
is detectable for A/B — 0.1 and A/B = 5.0, for exmaple, 
which are relevant to 1 H nuclei in MnCu(pbaOH)(H20).3 
with J/k B ~ 34 K [3] and 19 F nuclei in Mn(hfac) 2 NITiPr 
with J / k-Q ~ 460 K [4] , respectively. When we measure 
MnCu(pbaOH)(H 2 0) 3 and Mn(hfac) 2 NITiPr at 100 K 
and room temperature, respectively, the three-magnon 
relaxation is possibly dominant for H < 0.6 T in both 
compounds. 

There exist pioneering T\ measurements on the layered 
ferromagnet CrCl3 [24] and the coupled-chain antifcrro- 
magnet CsMnCi3-2H 2 [25], which give evidence of the 
relevant three-magnon relaxation process in accordance 



with the Pincus-Beeman formulation [18]. However, they 
are both low-temperature measurements under the ex- 
isting three-dimensional long-range order. A similarly- 
motivated experiment on the one-dimensional easy-plane 
ferromagnet CSN1F3 [26] was performed at higher tem- 
peratures but explained in terms of solitons, possibly be- 
cause of the lack of the spin- wave formulation in one di- 
mension. Here is a modern theoretical tool— the modified 
spin- wave theory based on the new scheme [12]. Let us 
observe the three-magnon nuclear-spin scattering in the 
one-dimensional quantum relaxation. 
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